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^ : ABSTRACT 

a^ : 

^ , We analyze the effect of pressure on the evolution of perturbations of an 

O , Einstein-de Sitter Universe in the matter dominated epoch assuming an ideal 

O ' gas equation of state. For the sake of simplicity the temperature is considered 

O ! uniform. The goal of the paper is to examine the validity of the linear ap- 

I . proximation. With this purpose the evolution equations are developed including 

{^X), quadratic terms in the derivatives of the metric perturbations and using coor- 

j> ! dinate conditions that, in the linear case, reduce to the longitudinal gauge. We 

^ . obtain the general solution, in the coordinate space, of the evolution equation for 

the scalar mode, and, in the case of spherical symmetry, we express this solution 
in terms of unidimensional integrals of the initial conditions: the initial values of 
the Newtonian potential and its first time derivative. We find that the contribu- 
tion of the initial first time derivative, which has been systematically forgotten, 
allows to form inhomogeneities similar to a cluster of galaxies starting with very 
small density contrast. Finally, we obtain the first non linear correction to the 
linearized solution due to the quadratic terms in the evolution equations. Here 
we find that a non null pressure plays a crucial role in constraining the non linear 
corrections. It is shown, by means of examples, that reasonable thermal velocities 
at the present epoch (not bigger than 10~^) make the ratio between the first non 
linear correction and the linear solution of the order of 10~^ for a galaxy cluster 
inhomogeneity. 



Subject headings: cosmology: theory, galaxies: clusters: general 
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1. Introduction 

The relativistic theory of the evolution of perturbations was initiated in 1946 by Lifshitz 
using a special coordinate condition known as the synchronous gauge. He linearized the 
Eintein's equations to obtain the evolution of perturbations. In fact, he found plane wave 
solutions for the radiation dominated epoch, assuming p = (l/3)p as equation of state, 
and for the matter dominated epoch, neglecting the effects of the pressure. The theory, 
with subsequent improvements, is referenced in many books of Cosmology (Peebles 1980; 
Zel'dovich & Novikov 1983; Landau & Lifshitz 1979; Weinberg 1972). 

However, the synchronous coordinate condition has two great drawbacks. The first one 
is consequence of the fact that it does not completely fix the coordinate system, allowing the 
existence of gauge modes. This problem can be handled using a gauge invariant version of the 
theory, started by Bardeen (1934) and collected by Mukhanov, Feldman & Branderberger 
(1992). This last review also shows an easier way to obtain gauge invariant equations using 
a coordinate condition that does not allow the existence of gauge modes. These coordinate 
conditions define what is known as the longitudinal gauge. 

The other inconvenience of the synchronous gauge is that the metric perturbation and 
the density contrast both depend on the second space-like derivatives of a potential. Then, 
great values for the density contrast imply great values for the metric perturbation, and 
in consequence the linear approximation in this gauge fails when the density contrast is 
bigger than unity. On the contrary, in the longitudinal gauge the metric perturbation is 
proportional to a potential while the density contrast is proportional to the laplacian of 
the same potential. So, the metric perturbation can be a very small quantity while the 
corresponding density contrast can achieve values greater than unity. For example, galaxy 
clusters develop a potential of the order of (f)/c^ < 10~^ varying at scales of i? ~ lMpc/6000/i; 
then, using the relation d = (1/6) A0 for the density contrast, in adimensional coordinates, 
we get 6 ~ l/(6i?^) ~ 10^ for h = 0.5. This makes, in principle, possible the validity of the 
linear approximation in the longitudinal gauge to study the formation of inhomogeneities 
similar to galaxy clusters. 

Then, the question arises why the linear approximation (linear in the metric perturba- 
tion) is always considered inaccurated to describe the evolution when the density contrast is 
bigger than unity. In this paper we are interested in analyzing when the linear approxima- 
tion begins to fail in describing the evolution of such objects. To do that, according to the 
previous paragraph, it is necessary to take into account here on that the spatial derivatives 
of the potential can be much bigger than the potential. Let us to point, in advance, that the 
pressure plays a crucial role in this issue. 
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So, in section 2 we write Einstein's equations in evolutive form, keeping quadratic terms 
in the first derivatives of the metric perturbation and neglecting quadratic terms in the 
potential. We use coordinates which simplify the evolution of the tensor components of 
the metric perturbation and which reduce, in the linear case, to the longitudinal gauge. 
Moreover, to complete the evolution equations wc need to give the stress tensor. In the 
matter dominated epoch and after the decoupling with radiation, the temperature of matter 
T{t) decreases as l/a^(t), where a{t) is the expansion factor, and the pressure becomes so 
small that it is usually neglected. But, as we shall see in this paper a non null pressure is 
necessary to keep valid the linear approximation. Then, we will consider the simple case of 
an ideal gas with an equation of state of the form p — {alTo/ma'^{t))p. Notice that, although 
the temperature is becoming very small, the evolution will increase the energy density and 
this is the reason for keeping the pressure. Under these conditions we write down, in section 
3, the evolution equations of the gravitational potential in the linear approximation. 

The evolution equations have two degrees of freedom: the potential and its first time 
derivative. In section 4 we find the general solution cf)^'^^ of the linear evolution equations for 
arbitrary initial conditions. Next, we use this result in section 5 to obtain the evolution of the 
density contrast and the macroscopic velocity starting from appropriated initial conditions to 
form an inhomogeneity similar to a galaxy cluster. The characteristic length of the structure 
will be given by the parameter e = r(l + ZiY^"^, as a sort of Jeans length, where Zi is the 
initial redshift and r — represents the present value of the random mean square 

velocity (r.m.s. velocity). 

Finally, we face the problem of the validity of the linear approximation. The validity 
criterion, working in the longitudinal gauge, cannot be based on the value of the density 
contrast, as we have commented above. Instead, it should be based on the value of the 
non linear corrections of the evolution equations. So, in section 6 we estimate the first non 
linear correction, (^^^ ^ due to the quadratic terms in the Einstein's equations, and obtain 
an upper bound estimation for the quotient F = \(\)^^^ j (\)^^^\. The linear approximation will 
be considered suitable if this quotient is small, although the density contrast has reached a 
great value. With this validity criterion, an inhomogeneity similar to a cluster of galaxies 
could be described with the linear approximation because we obtain a F of the order of 10~^. 
Notice that, if we take the well known solution for p = 0, we obtain F bigger than unity. So, 
the pressure plays a crucial role for the validity of the linear approximation. 
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2. The evolution equations 

We assume in this paper the concepts and notations usual of the 3 + 1 formaUsm of 
general relativity (Smarr & York 1978). We shall consider a perturbation of an Einstein-de 
Sitter Universe, so we put the metric in the form: 

ds^ = -a^{(t))dt^ + -fijdx'dx^ (1) 

The shift vector (3^ has been taken null, and the lapse function a will be choosen conveniently 
later. We write the tridimensional metric in terms of a scalar and a trace-less tensor: 

-fij ^ a\{l - 2(l))Sij + aij) (2) 

where a — a{t) denotes the scale factor of the Einstein-de Sitter Universe, and aij is a 
tridimensional tensor verifying S'^^'amn — 0. In the following we shall neglect quadratic 
terms in the metric perturbations, (f) and 0"^, but those which are quadratic in its first 
derivatives (as was explained in the introduction). We shall need the extrinsic curvature of 
the surfaces t ^constant, 

1 1 

Kij := -—dt-fij = {{H{1 - 24>) - dt4>)5ij + Haij + -dtaij], 

la a 1 

being H — a/a the Hubble constant, and the Ricci tensor of the tridimensional metric ^if 
Rij = (1 + 20)0,,,- + 30,,0,,- + ((1 + 20) A0 + (V0)2)5,,- - ^Aaij - 5^^a^im,mj) 

where the operators A and V are referred to the euclidean tridimensional metric. 

Splitting the energy tensor in parallel and orthogonal components to the vector field 
u = {l/a)dt, 

T^^ = puV + ph"" + qf'u'' + q^u" + n"", (3) 

one gets the corresponding energy density, flux of energy and stress tensor. Then, a Cauchy 
problem with constraints can be stated in General Relativity (Bruhat & York 1980). Over 
a space-like surface t = ti a tridimensional metric and an extrinsic curvature tensor are 
suposed to be given. These tensors evolve in time according to the following equations: 

dtKij = -DiDja + a{Rij + trKKij-2KiaK^)+4nGa{p- p)-fij-8nGanij ^ ' 

with trK representing the tridimensional metric trace of the extrinsic curvature tensor and 
Di the tridimensional covariant derivative. In these equations the components of the stress 
tensor, p and tt^-, must be choosen from the beginning, as we will do below. 
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The energy density and the flux of energy are linked by constraint conditions to the 
Ricci tensor and the extrinsic curvature: 

IGTrGp = (trKf - tr{K^) + R 

SnGqi = D-Kai - DitrK ^ ^ 

where R is the scalar curvature of the tridimensional metric. If one knows at some ini- 
tial instant ti the flux of energy and the energy density, one must solve first the constraint 

equations (5) to determine a valid set of initial conditions ^j*j{x),K*j{x). Then, the evolu- 
tion equations (4) determine the '~fij{t,x), Kij{t,x) for t > ti. Substituting them into the 
constraint equations one gets the evolution of the energy density and the flux of energy. 

Finally, we shall assume a one-component Universe with a pressure tensor of the form: 

= + 7r(f(0,t) (6) 

p = pB + EA(t) + p^^\(j),t) (7) 

where means trace-less component, A and E are only functions of time and p^'^^ and 

(2) ' 

are a scalar and a second order 3-tcnsor formed with the 3-vectors , crim,m and its first 
time derivatives respectively. This assumption is quite general because it allows to consider 
an ideal gas as well as solutions of an Einstein- Vlasov problem. Next, we shall develope the 
evolution equations taking into account these last expressions. 

Let us start splitting the second evolution equation into trace-less and trace part equa- 
tions. The trace-less part is: 

where 

Sij = + 20) -a' - S'KGaA)[(t)^ij]^ + {-a" - 2a' + 3a)[(/),i0j]* - STrGaTrJ^^ (9) 

An appropriate election of the lapse function q;(0) can simplify the problem. Lifshitz used 
the Gaussian gauge: a — 1, but this choice has more than one inconvenience, as we have 
pointed out at the introduction. Looking at equation (9) we see that the best choice is to 
take a such that the coefficient of the Hessian vanishes. That means to take a as the solution 
of the equation: 

a(l -I- 20) - a - S^GaA = 

which is « = e^i'^^'^^, with bi = 1 — SnGA. With this election, the coefficient of [0,i0j]* in 
the expression of S^j becomes —2 -|- 327tGA, and the evolution of the trace-less component 
results: 

(Jj (Jj (Jj 
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As to the trace component of the evolution equation, it writes down as: 
dl<j> + AHdtcj> - (^1^ + 47rG£;)A0 + ^d.drr^a.m = \{dtct>f - ^(V0)' + A^^Gp^^^ 

We have found convenient to introduce the conformal time 7], defined by dt = a^drj. In 
this time coordinate the expansion factor writes down as a^rj) = a^rf with related to the 
Hubble constant by Oo = 2/i7o. Then, the final form of the evolution equations is: 

dl<i> + %^<i> - ^(2A + 3Ea2)A0 + ^a^„ = \{d^(t)f - ^{V<i>f + AnGa^p^^^ (10) 

d^aij + Ur^a,, - Aaij - 2[a(,>„,,„)]* = 4(1 - IGTrGA) + 167rG(l - 87rGA)7rlf (11) 
In addition to these equations, the constraint conditions should be considered: 
1 3H 3H^ 5 3 1 

AnGp = -A0 - —0,^ + - 2^1'^) + ^C^*^)' + + ^^irn,im (12) 

AnGqi = -hiH(j)^i - -0,^^ + ^aim^mri + Haim,m (13) 

So, it remains to give the functions A, E, p(^) and tt^J'' appearing in the pressure tensor. 
As we are interested in the matter dominated epoch, it is reasonable to consider the one- 
component Universe as an ideal gas, with energy tensor T^^, — pcW^Wy + Pc{g^lu + w^Wu), 
equation of state Pc — {T/m)pc and temperature evolving as T = const /c?. In the Appendix 
A we show how this assumption means to take ^4 = and E Ri r'^/AnGa'^, where r = y^T^Jm 
is the r.m.s. veloctiy at the present epoch, and neglect the second order expressions p^^^ and 



(2) 
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Let us repoduce the complete set of equations: 

a^c/, + - ^ A0 + la,^,,^ = l(a^0)2-i(V0)2 (i4) 

4 

d'^aij + -dr^aij - Aaij - 2[a(jm,im)Y = (15) 

5 = —Act) - Tl(t>,'n - 20 + —(Tim,im (16) 



A-nGqi = -Hcp^i - -0,^^ + ^aim,mv + Haim,m (17) 



where we have substituted the energy density p by the density contrast 5 using the relation 
5 — {p — Pg)/ Pg, with Pg — SH'^/SnG the background energy density. 
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2.1. The longitudinal gauge 

The longitudinal gauge, unlike the Gaussian gauge, fixes definitely the coordinates. 
This is a well known fact, but let us give here an argument, which may be useful for other 
purposes. 

If we start with the Robertson- Walker (R-W) metric in canonic coordinates, ds^ ~ 
—dp + a?{t)5ijdx'^dx^ and we introduce new coordinates: 

t — t + ip{t,x) 
x' = x^ + ^i{t,x) 

and impose the coordinate conditions goi — 0, one puts the metric in the form: 

ds'^ = -(1 + 2ip)dP + a^{t)[{l + 2H^p)5mndx^dx^ + 2iirn,n)dx'^dx''] (18) 

with 

1 [I 

Cm{t,x) = -tp^rn + Cm{x), = 2 / -i:^dt 

2 J 

being ip the time derivative of and Cm arbitrary functions of the space-like coordinates. 
Then, the lapse function and the tridimensional metric are: 

a = 1 + 



iT - At) 



(1 + 2H^ + + \Cm,m)5ij + + KUn) 



Comparing these expressions with equation (2) we obtain the gravitational potential and the 
trace-less tensor as: 

Consequently, if we choose the longitudinal gauge, i.e. a — 1 + (f), the function (p, which 
define the new time coordinate, should satisfy the equation: 

d^<p + ^A<p + H<p + Hdt<p = 0. (19) 

Notice that the coefficients of the second time derivatives and of the laplacian operator in 
this equation have the same sign. This makes impossible to construct a time-like foliation 
with the function (p, apart from the case (p — 0. So, we can conclude that if we develope a 
R-W perturbation in the longitudinal gauge, it is impossible to recover a R-W space-time in 
other coordinates. 

However, if one chooses the Gaussian gauge, i.e. a — 1, the function (p should now 
satisfy the equation (p — 0, making possible to introduce new time coordinates. This forces 
to characterize the gauge modes, i.e. the false metric perturbations. 
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3. The linear evolution equations 

We shall study now the linear equations: 

6 1 

^^(|)+-^r,(|)-^A(j)+—a,m,^m = (20) 

4 

d^aij + -dr,aij - Aa^ - 2[a(jrn,im)Y = (21) 

2 2 

S = |-A0 - - 20 + (22) 
AirGqi = -H(f)^i 0,^j + —aim,mri + Haim,m (23) 

Firstly, we consider that the trace-less symmetric tensor cjjj can be decomposed (York 1973) 
in a transverse part cr^, verifying aj^ ,^ = 0, and a longitudinal part afj, with aj^^ „^ = aim,m- 
These two components are orthogonal and evolve independently. Then, one can distinguish 
three modes in our problem: the scalar one 0, the transverse tensor aj- and the longitudi- 
nal tensor 0"^^. The constraint equations show that the scalar mode is the most important 
contribution to the density contrast (22) and to a rotational-free flux of matter (23); the 
longitudinal tensor mode contributes weakly to the density contrast and, what is more inter- 
esting, it is the only possibility of producing a non null rotational component of the velocity 
field. Relative to the tensor transverse mode, it does not contribute nor to the density neither 
to the energy flux, in fact it represents the emission of gravitational waves. 

Next, given that the scalar and the longitudinal modes are coupled and that the double 
divergence of the longitudinal part appears in the evolution of 0, we can rewrite the equations 
introducing the scalar 9 = — (l/12)(TiTO_j^. In this manner, the evolution of gives the 
evolution of 9 and the weak coupling between the scalar mode and the longitudinal one 9 
is given by: 

+ -a,0 - 4 A0 = 9 (24) 

dl9 + -dr,9 -Ia9 = (25) 
' 77 3 

Both equations arc hyperbolic, but there exist a great difference between them due 
to the time dependent coefficient rj'^ that appears in the evolution of (p. This coefficient 
causes that the characteristic curves of the first equation do not scape to inflnity as do the 
charateristics of 9, because in this case the laplacian operator has a constant coefficient. This 
can be seen clearly in the case of spherical symmetry, where the charateristic curves r{'r)) of 
(24) tends to a finit limit as the conformal time tends to infinity; while for the variable 9, the 
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characteristic curves of (25) scape to infinity. This makes possible that a hnear hyperbohc 
equation might describe the increasing of density in bounded regions. 

Morever, this difference makes the couphng between the scalar and longitudinal modes 
almost irrelevant, because small initial values for ^ in a small region disperse to infinity. 

Finally, we observe an important difference between our evolution equations an those 
in the Gaussian gauge (lapse function a — 1 and shift vector (3 — 0) used by Lifshitz. In our 
gauge, the evolution equation for the scalar mode can be reduced to a unique equation for 
a unique function, while in the Gaussian gauge the scalar mode is described by two coupled 
equations. 

In the next section we shall find the general solution in the coordinate space of the 
evolution equation for the scalar mode, neglecting the coupling with 9 or assuming ^(?7j, x) — 
dr,6(r]i,x) = 0. 



4. The Cauchy problem of the linear evolution equations 

From equation (24) and in the case of null initial conditions for 9, we can consider the 
following initial value problem: 

d',<l> + ljdJ-^A<l> = ^26) 
HVi, = M^) , dr,(l){r]i, = (l)'i{x) 

with (pi{x) and 0-(a:) two arbitrary functions. We shall solve this problem using the method 
of Fourier transforms. 

Firstly, let us denote by (j){r].s), with s G M^, the Fourier transform of (j){r],x) with 
respect to the spatial coordinates. In the Fourier space, equation (26) transforms into an 
initial value problem for an ordinary differential equation: 

dl^ + ^d,^ + ■ = ^2^^ 

(j){rji, s) = ^i{s) dr^(j){rji, s) = (j)-{s) 

So, the first task is that of constructing a system of fundamental solutions, which consists 
on two solutions (f)i{r], s), (t)2{'rj, s) that satisfy the initial conditions: 

0i(?7i,s) = 1 , <9^0i(?7i,s) = 
hiVi, s) = , dr,(i)2{rii, s) = 1 
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These fundamental solutions can be obtained using complex Laplace transforms (Smirnov), 
having the following result: 

^ 3 f sin kg cos kg\ r)i(3ri-r)i) sin kg nf 

MV, s) = -, - g^ j + ^ + ^ cos kg (28) 



1 / \ 9r]i f sin kq cos kq osin kq \ , 

MV:S) = -ki^-^-g^-g^^) + 



(29) 



37]- f sin kg _ cos kg \ rjf sin kg 

^ ^yn^ ^n^)^^ k 

where k stands for the modulus ol k — y/s ■ s, e — r/rji, and g = e{l — ^). Then, the 
solution of (27), in the Fourier space, is expressed in terms of the fundamental system as: 



4>{V, s) = $i{s)4>i{r], s) + ^l{s)4>2{r], s). (30) 

The next task is to obtain the Green's functions Qi{r],x) and Q2{fli^)i defined as the 
inverse Fourier transform of the fundamental system {0i(r;, s), 02(^7,-5)}. As we show in 
Appendix B, these Green functions are: 

= ^.(,-,.(5^-^)M^.4,a,(M^) (31) 

- (i;^(^^-^vj|,)«(.-o+4.^ (3.) 

Therefore, the solution of (26) in the coordinate space is expressed as the convolution product 
of the Green's functions with the initial conditions: 

(j){r], x) = Qi{r], x) * 0j(a;) + Q2{r], x) * (^[{x) (33) 

where * stands for the convolution product with respect to the spatial coordinates. Looking 
at the Green's functions we can observe that the solution tends rapidly to a limit when the 
conformal time tends to infinity, this limit has a simple expression: 

Hoc, x)^-^f 0,(e)cie + ^ / (e'- 1 ^ - e mm^ (34) 

as tridimensional integrals of the initial conditions. The e parameter in expression (34) can 
be also written as e = T^/T+~Zi, with r the r.m.s. velocity of the matter component at the 
present epoch and Zi the initial redshift (recall the relation 1 + z = 1/77^ in an Einstein-de 
Sitter Universe) . This parameter will be crucial to the study of evolution because it will fix 
the characteristic length of the evolved structures as a sort of Jeans length. 

In the case of spherical symmetry the convolutions reduce to unidimensional integrals, 
whose expressions are obtained in the Appendix B. They will be used in the sequel to discuss 
some examples. 
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5. The evolution of fluctuations with non nuU thermal motions 

Having got the general solution of the linear initial value problem, the constraint equa- 
tions (22) and (23) determine the density contrast and the flux of matter as simple functionals 
of the metric perturbations. In the linear approximation we have that these equations for S 
and Qm reduce to: 

6{rj,x) = ^A(/)-20-#' (35) 

A-KGaqm = -aH(f)^rn - 0,r,m (36) 

where we only have to substitute the expression (33) of the solution 0. Notice that in the 
case of an ideal gas, the flux of energy in the longitudinal gauge represents the macroscopic 
velocity: Qm = pVm- In particular we are interested in the velocity norm, whose expression 
is: 

l^(^'^)l7=I^|V(|0+^0,.)|- (37) 

Let us study these expressions when the general solution (33) is substituted. We shall 
assume that at some initial redshift Zi we know the initial conditions for the potential (t)i{x) 
and its first time derivative (l)'i{x). We shall begin with a qualitative description of the 
evolution of the density contrast based on the following reduced expression for 5: 

S{n,x) « '^iQi{oo,x)*A(Pi{x)+Q2{oo,x)*<j>[{x)) (38) 



Qi{oo,x) = -l-^H{e-r) (39) 

Q2{oo,x) = ^(e'-r')H{e-r) (40) 

where we have only considered the laplacian term neglecting the (f) and 0' contributions and 
we have also taken the asymptotic values for the Green's functions. The idea is to obtain 
estimations of the convolutions using the Holder inequahties (Hormander 1989). In this case 
we have enough with the relations \\ f * g ||oo^|| / ||i|| 9 Woo, where || / ||oo means sup \ f |, 
and II / 111 means / | / | dx. We obtain in this way, ior r) >> r)i, two upper bounds for S: 



2 

S{r},x) ||oo< ^( II Qi{oo,x) lloo II A0i(x) 111 + II (52(00, x) lloo II 111 ) (41) 



II S{r},x) ||oo< II Qi{oo,x) 111 II A(f)i{x) \\^ + \\ Q2{oo,x) ||i || A(f)'-{x) \\^ ^ (42) 

This will allow to reach the main conclusions with little calculations, treating separately each 
one of our two degrees of freedom (f)i{x) and 0^(x). In the subsections we shall give more 
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details in numerical examples where we shall assume spherical symmetry. To do it we need 
to fix the parameters of the problem. We have three parameters related by the condition 
e = T\J\ + Zi. As we have mentioned above the final characteristic length of a structure will 
depend on the value of e, so if we want to discuss galaxy clusters it will be convenient to take 
e ~ IMpc/aoi tto = GOOOh^^ Mpc, which is of the order of an Abel's radius for this value of 
e. And assuming Zi = 5000 we obtain r ~ 10^^ for the r.m.s. velocity at the present epoch. 
Notice that this value can be supported by hot particles as neutrinos with non null mass. 

In these examples, to identify the final structure as something similar to a galaxy cluster, 
we shall estimate the total mass at the present epoch r] = 1 and inside an Abel's radius 
Ta = 1.5h^^Mpc as a function of the amplitude A and the initial characteristic length R: 

M{< r„) = 47r2.7 x IO^^Mq f " S{r, 1, R, A)dr (43) 

Jo 

Recall that a typical galaxy cluster has M(< r^) ^ 3h~^ x 

In the subsections we examine the possibility of generating structures similar to a galaxy 
cluster starting from reasonable initial conditions for the potential, that is, such that the 
initial density contrast and the macroscopic velocity be small. We shall also require that 
in both examples the gravitational potential at the decoupling of matter and radiation be 
smaller than 6 x 10~^. 



5.1. Initial conditions of the form (pi 0, (f)'- — 



Let us start studying the case of 0i 7^ and 0' = 0. In this case, from equations 
(35) and (37), one gets the initial density contrast and the initial macroscopic velocity as: 
Si{x) ~ ^A0j(a;) and | V{rii,x) y | V0i |. Let us consider an initial potential of the 
form (f)i{x) — — A{1 + r'^ / R'^)~^^'^ , in which we have two parameters, the amplitude A and the 
characteristic length R. The laplacian of this function is A0j = 6S(1 + r^/i?^)~^/^, where B 
is given by the relation A — 2BR? and is proportional to the initial central density contrast, 
5i(0) = rifB. Under these considerations, from equations (41) and (42) we have that the 
final density contrast is upper bounded by: 



5{r],x) \\ao< minor of 



So, one distisguishes two cases: 



If < e we have || S{rj,x) \\< (?7/77j)^(i?/e)^5j(0). In this case the density contrast 
decreases until the instant 77*, for which {r]^/r]if{R/ef — 1, and then begin to increase 
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with a power law. This possibility was stated first by Gilbert (Gilbert 1966) using a 
Newtonian approximation to the Einstein equations. 

• If i? > e the object grows from the very beginning, tending to a power law: || S{r], x) \\< 

Getting on with the example we will assign numerical values to the free parameters 
A and R. We shall choose them such that the total mass inside an Abel sphere of radius 
Ta = 1.5h~^Mpc be about 3h~^ x 10^^ solar masses and, at the same time, keeping bounded 
the gravitational potential by | (f){ri,x) |< 6 x 10~^. For example, taking A = 2.8 x 10~^ we 
can obtain the mass M{< ra) as a function of the characteristic length R, whose graph is 
given in Figure 1. 



M(<ra > 




Fig. 1. — The mass M(< r^) inside an Abel's radius (r^ = 1.5h~^Mpc ) at the preset epoch 
as a function of the initial characteristic length R, given by equation (43) when we evolve the 
density contrast taking into account only the first degree of freedom, i.e., having as initial 
conditions 0i = -^(1 + r^/R^)~^^'^ and 0- = 0, with A = 2.8 x 10~^. We assume the values 
e = IMpc/ao , Oo = 6000h~^ Mpc , h — 0.5 and Zi — 5000 for the e-parameter and the initial 
redshift. The r.m.s. velocity at the present epoch parameter takes the value r = 1.2 x 10~^. 

This figure shows that we can choose R — 1.6e to obtain the mass of a typical galaxy cluster. 
With this election, equations (37) and (38) give the evolution of the velocity and the density 
contrast, whose graphs at the initial and final time are shown in Figure 2. 
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Fig. 2. — Initial (left figures) and present (right figures) values for the density contrast and 
the macroscopic velocity, evolving the first degree of freedom under the same conditions as 
in Figure 1 and with R — 1.6e. The mass inside an Abel's radius Va — 1.5h~^Mpc at the 
present epoch is about 6 x IO^^Mq. 

In Figure 3 we show the evolution of the central value of the density contrast in the 
case where R < eto illustrate the bouncing of inhomogeneities with typical length R smaller 
than e, as we have described above. 

As conclusion of this subsection we can say that the evolution of Einstein's equations 
with initial conditions of the form {0i(a:) 7^ O,0-(x) = 0} is equivalent to the evolution 
goberned by Newton's equations (Gilbert 1966) with initial conditions {6i{x),Vi{x)} given 
at the beginning of the subsection. The second example will lead us to a quite different 
conclusion. 
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Fig. 3. — The central density contrast obtained evolving the first degree of freedom. The 
right figure shows the evolution under the same conditions as in Figure 1 and with R— 1.6e. 
The left figure shows the bouncing produced on the central density contrast when we take 
an initial characteristic length R — 0.2e. Then we can see that the e parameter plays the 
role of a Jeans length. 



5.2. Initial conditions of the form 



0, 0; 7^ 



Next we are going to study the case where (pi = and 0^ 7^ 0. In this case, the initial 
density contrast and the macroscopic velocity are not given by the initial potential but by 

2 

its initial first time derivative in the form: 5i{x) — —r]i(f)[ and | V{r]i,x) y | |. Now, 
wc shall consider an initial condition of the form = —A{1 + r^/i?^)~^/^, with A and 

R two free parameters. Defining A — 2BR^ we can write A(/)^ = 6S(1 + r'^jR?)'^!'^ having 
that the final density contrast is upper bounded by: 



^{r\^x) \oo< minor of 



1 o R^ 1 2 ' 

-7] Vi-^B , -7] TJiB 



Let us remark the main difference with the previous case. If we consider here R > e, 
the final density contrast will be bounded by {l/5)riiB, but now B is not constrained to be 
small because the initial density contrast and the macroscopic velocity as functions of B and 
R are given by: 

I V{Vi,r) \j 

having that for objects much smaller than the horizon at the present epoch {R = IMpc/ag 
with Oo = 6000h~^ Mpc imphes R ~ lO"*^), rjiBR^ can be small even when B » 1. There- 



2n,BR^ 
(l + rVi?2)V2 

rjfBr 
3(1 + rVi?2)3/2 



(44) 
(45) 
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fore, the second degree of freedom allows to reach great values of the density contrast starting 
with very small initial density contrast. 

As in the previous subsection, we have to take values for A and R such that the total mass 
inside an Abel sphere of radius Va — 1.5h~^Mpc be of the order of 3h~^ x 10^^ solar masses 
and, at the same time, keeping bounded the gravitational potential by | (f){r], x) |< 6 x 10"^. 
For example, taking A = 3 x 10"^ we can determine M(< Tq) as a function of R, whose 
graph is represented in Figure 4. 

MC<rJ 




Fig. 4. — The mass M(< Va) inside an Abel's radius (r^ = 1.5h~^Mpc ) at the preset epoch 
as a function of the initial characteristic length R, given by equation (43) when we evolve 
density contrast taking into account only the second degree of freedom, i.e., having as initial 
conditions 0i = and 0^ = —^(1 + r^/i?^)"-*^/^, with A — 10"^. We assume the same values 
for e, T and h as in Figure 1. 

From this picture we obtain that a good value for the characteristic length is R = 1.6e. 
These values of the parameters allow us to evolve the density contrast and the velocity, 
whose evolution is represented in Figure 5. 

Then, we can conclude that the second degree of freedom, taking an appropriate value 
for 0', allows the formation of great structures starting from very small initial values for 
the density contrast. This case has no Newtonian analogue because now we have significant 
initial time derivatives of the gravitational potential. In other words, a Newtonian evolution 
starting with initial density contrast and macroscopic velocities as given by expressions (44) 



-17- 



■ IQQQ Vj. ■ IQQQ 




Fig. 5. — Initial (left figures) and present (right figures) values for the density contrast and 
the macroscopic velocity, evolving the second degree of freedom with the choice of parameters 
done in Figure 4 and with R — 1.6e. The mass inside an Abel's radius = 1.5h~^Mpc at 
the present epoch is 6 x IO^^Mq. 

and (45) will produce a structure with a very small final density contrast at the present 
epoch. 

Finally, let us remark that the second degree of freedom can be described geometrically 
as follows: the initial surface rj — rji has null intrinsic curvature (null laplacian of 0i) and 
highly inhomogeneous extrinsic curvature (great space derivatives of 09 . 

6. On the validity of the linear approximation 

In this section we come back to the non linear equations (14) and (15) in order to study 
the validity of the linear approximation. Introducing the function 9 as in section 3, we obtain 
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a coupled system of evolution equations for the couple of functions (0, 9): 

+ -a,(/> - 4 A0 = ^ + l(9,</,)2-l(V0f 
' r) r)^ 2 b 

with initial conditions (f){rji,x) = (f)i{x), 9^0(r/i,a;) = 4>i{x), 9{r)i,0) = 0, and dr^9{r)i,x) — 
0. This is a semilinear hyperbolic initial value problem. In the Courant-Hilbert's book 
(Courant & Hilbert 1962) the unicity of solutions of this kind of problem is shown by means 
of the convergence of iterations. This supports the fact of considering the first iteration 
as criterion for validity of the linear approximation. In the following we shall focus on the 
reduced equation: 

+ -9,0 - -,A<t> = -liW<py (46) 
' r) T)^ b 

because it contains the essentials of the problem. Let us denote by (p^'^^ the solution (33) to 
the linearized equation and by (p^^^ the first non linear correction, namely the solution of 

9^0(1) + _ I!a0(i) = -i(V0(°))^ (47) 

77 77* 6 

with null initial conditions. Using the Fourier transform method, the problem reduces to an 
ordinary differential equation: 

+ + ■ = (77, k) (48) 

where L^^\r], k) stands for the Fourier transform of the quadratic term — |(V0*^°-')^, which is 
easily solved by the method of constants variation. The fundamental solutions {01(77, s), 02(^7, -s)} 
of the homogeneous equation were obtained in section 4, see expressions (28) and (29). Then, 
the solution of (48) can be expressed by means of integrals: 

and in the coordinate space it results: 
0^(7^,2;) = f 



(r9{ri) 
Q^*J^ il-irQi*L(^H9,x)dg- 

(l-f)U*L(0)(<7,x)dJ 



(49) 
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being Qi and Q2 the Green's functions given by (31) and (32). To decide about the vahdity 
of the hnear approximation we need to compare | (p^^\rj,x) \ with | cf)^^"' (7], x) \ . This is not 
easy to do directly, but using norms we shall get an upper bound for the first non linear 
correction which will be enough to discuss the problem. Then, we have: 



|0«(r7,x)| < !Ill^\Q,* J^'^'\l-l)^Q,:,L^^)(g,x)dg\+... 

where dots means the same expression but interchanging Q2 for Qi- As before we are 
going to use the Holder inequalities, in particular ||/*5'||oo < ll/lli lls'lloo ^-nd 
II / * 5' II2 < II / 111 II 9 II2, where || g II2 means / | 5' |^ dx. We get in this way: 

U^'Hoo,x) lloo < ^(^\\Q2{oo,x) \U J\l-^f\\Qi{g,x) \U\L^''\g,x) ||ici^ + ...) 

In the sequel we shall obtain this upper bound for the numerical example studied in 
subsection 5.2 corresponding to the initial conditions of the form (f)i{x) = and 0- 7^ 0. 
So, we write (j)^^\ri,x) = Q2{r],x) * (^^{x) and get the estimation || L^'^\g,x) ||j < ^ || 
Q2{g,x) 11^ ^a=i II Va0^(a;) H^. Substituting this into the previous equation we obtain: 

1 ^ 

II 0«(oo,a;) lU < -G{r]„e)J2\\ ^a<PK^) \\l (50) 

a=l 

with 

Givi,e) = j(^\\Q2{oo,x)\U \\Q,\\^ \\Q2\\ldg+ 

+ II gi(oo,a;) 111 ^ (1 - f)' II Q2 lloo II Q2 111 dg^ 
Evaluating the corresponding LP norms we have: 

which reduces, once calculated the integrals, to: 

G{m.e) = l^-^'^, (51) 

Then, in the case considered in the previous subsection 5.2, a simple calculation gives 
Sa=i II ^a^i(^) II 2 — 37r^S^i?^. Substituting these results into equation (50), we get that 
the upper bound C^^^ to the first non linear correction is: 

I (t>^^\oo,x) I < II (t>^^\oo,x) lloo < C^^^ = 10-V^S^i?2 (52) 



e3 
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Now, we have to compare this bound with the norm of the hnear solution (/)(°). To 
do it we form the quotient F = C^^V | 0^°^(oo,O) |. Given that for this example we have 
spherical symmetry, we can use the unidimensional integrals of Appendix B to calculate the 
modulus of . Then, T expresses as a function of the initial characteristic length R, which 
is represented in Figure 6. 

r 



Q,QQQ8 
Q,QQQ6 
Q ,QQQ4 



Fig. 6. — Validity of the linear aproximation of the evolution of the second degree of freedom 
showed in Figure 5. The function F is an upperbound, at the present epoch, of the ratio of 
the first non linear correction and the linear solution for the gravitational potential. 

As we can see in this figure, the quotient F is below 10~^ in the range R < 2e, then the 
relation | 0*^^)(oo,x) |< 0.01 | (j)^^\oo,x) \ is verified. So, as we have choosen R — 1.6e, we 
can neglect the first non linear correction and consequently, the linear approximation is an 
accurate description of this problem even having reached a great final density contrast. 

Therefore, we can conclude that a thermal velocity r of the order of 10"^ at the present 
epoch makes possible to follow with the linear approximation the formation of an inhomo- 
gencity similar to a galaxy cluster. We have seen also that the linear approximation comes 
into problems with smaller values for the thermal velocity. 

7. Conclusions 

In the current theory on evolution of perturbations, the matter dominated epoch is 
considered as a fluid with null pressure, and the evolution is described using the linear 




- 21 - 



approximation until the density contrast becomes of the order of unity. However, as we 
have shown in this paper, if pressure is properly considered, the evolution with the linear 
approximation can be extended to values of the density contrast bigger than unity. We 
have assumed an isotropic pressure of the form p = {a'^To/'ma'^{t))p, which corresponds to an 
ideal gas with uniform temperature, and we have used reasonable values for the temperature. 
Concretely, in our examples we have taken a random mean square velocity of the order of 

In the following we summarize the main steps we have followed to get these conclusions: 



1. We have stated a Cauchy problem using sistematically the 3+1 formalisme of General 
Relativity and neglecting quadratic terms in the metric perturbation. But, given that 
for inhomogeneities at scales of a few Mpc the spatial derivatives of the potential 
are much bigger than the potential, we have kept the quadratic terms in its first 
derivatives. These non linear corrections will only be used to study the validity of the 
linear approximation (see point 3). As usual, the coordinates are fixed by choosing 
the lapse function a and the shift vector (3. We have put /3 = 0, and a = e'^^^^'^^ has 
been taken in order to simplify the evolution equations. When linearized, this choice 
of coordinates is called the longitudinal gauge. 

2. We have obtained the solution of the linearized Cauchy problem for a one-component 

Universe in the matter dominated epoch and assuming an ideal gas equation of state 
p = {alTo/ma^{t))p. We have expressed this sohition in terms of convolution integrals 
of the initial conditions, which in our case arc the initial potential (j)i{x) and its first 
time derivative 0^(3;). We have also studied how to obtain a density inhomogcncity 
similar to a cluster of galaxies, i.e., how to get a total mass inside an Abel radius of 
the order of 3/i~^ x 10^^ solar masses. We have considered separately both degrees of 
freedom obtaining in both cases an inhomogeneity similar to a galaxy cluster at the 
present epoch. But the initial density contrast and the initial macroscopic velocity in 
each case are very different: 

(a) With initial conditions of the type {(t)i{x) ^ , (t)'i{x) = 0}, see subsection 5.1, one 
gets a cluster of galaxies starting at redshift Zi — 5000. These inital conditions 
correspond to an initial density contrast of Si ~ 0.1 (which becomes about 0.5 at 
the recombination redshift), and a macroscopic velocity of | |~ 3 x 10~^. The 
results of this case can be also obtained with a Newtonian analysis starting with 
the same density contrast and macroscopic velocity. 

(b) With initial conditions of the type {0i(x) — , 0-(x) 7^ 0}, see subsectio 5.2, one 
gets a galaxy cluster starting at the same redshift Zi — 5000. But now one has a 
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very small initial density contrast 5i ~ 0.0001, and a similar macroscopic velocity, 
I Vi |~ 8 X 10~^. Unlike the previous case, this evolution has no Newtonian 
analogue. 

The second degree of freedom, which is currently forgotten, may rapidly produce inho- 
mogeneities similar to galaxy clusters starting from faint initial density contrast. On 
the contrary, the first degree of freedom needs an excessive initial density contrast. 

3. We have estimated the first non linear correction to the linear approximation and used 
it as the criterion of validity. Our first results seem quite interesting: assuming at 
the present epoch a thermal velocity of the order of r 10~^, the quotient between 
the first correction to the gravitational potential 0^^^ and the linear solution is 
small than 10~^. Therefore, we can conclude that the linear approximation is an 
accurate description of the formation of a big structure if the effect of the pressure is 
not neglected. 



A. Description of an ideal gas in the longitudinal gauge 

An ideal gas in the matter dominated epoch with four-velocity w has a perfect fluid 

energy tensor T^i, = pcW^Wy + pdo^iu + w^Wu), with equation of state pc = Tpc/m being 
T = T^a^^ja^^ where is the temperature at the present epoch. The four- velocity w is 
related to the four- velocity u = {l/a)dt hj w = '~f{V){u + V)) where V is the macroscopic 
velocity of the Einstein-de Sitter perturbation in the longitudinal gauge, which is given by: 

47rGp a 2a 

where p is the energy density for the observer u. Using this relation, the energy tensor 
transforms into T^^ = pu^u^ + p{g^^ + m^m^) + q^jj,^ + g^^M^ + U^^, where now p, p and U^^ 
are quantities referred to the observer u and are given by: 

p^p, + 0{V') (Al) 
p^Pc + ^pV^ + OipcV^) (A2) 



a 



n,j=p[ViVj]' + 0{pV^) (A3) 

Given the equation of state and the relation pc — Pg{l + S) we can write Pc — Tpg/m-\- 
5pgT/m, which allows us to get: 

T / 1 2>H \ 
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where we have also introduced the r.m.s. velocity at the present epoch = {To/m). Prom 
this expression and comparing with equation (7), we identify: 

As for the anisotropic pressures, in the same way, equation (A3) gives: 

and comparing with (6) we obtain A = 0, and irf^ — Hij. 

Let us write the non linear evolution equations (10) and (11) as follows: 

C\^,a) = l(a,0)2-i(V0)2 + 47rGaV2) 



C\a) = 4(1 - 167rG'A)[0,i0j]* + 167rG'(l - %t:GA)t:\j 
where we have used a compact notation for the first members of the equations. These 

^3 



expressions, substituting A, p^"^^ and 7r|^'', transform in: 



L\4>, a) = l{d,<pr - 1(V0)2 + ^^^^^ (H\V<pr + 1(V0')^ + 2^V0 ■ V0') (A7) 
C'ia) = 4[</,,0,,]* + ^^^^^ ^H'[<l>,<f>jY + + + M',Y^ (A8) 

And taking into account that the density contrast S will have a great value in the structures 
we are interested on, the terms where it appears can be neglected. 

Let us to point out that keeping these terms only would produce small corrections to 
the upper bounds estimations for the first nonlinear correction to the linear approximation 
as obtained in section 6. 



B. Obtaining the Green functions 



In this appendix we are going to summarize the process to obtain the Green's functions 
associated to the general solution of our initial value Cauchy problem (26). This general 
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solution in the Fourier space has the form: 

4>{V, s) = 4)i{s)4)i{ri, s) + (f)[{s)(f)2{ri, s). 

being and (pi{s) the Fourier transform of the initial conditions and (f)i{r], s) and 02(?7, s) 
the fundamental solutions given by (28) and (29). In order to have the general solution in 
the coordinate space we need to make the corresponding inverse Fourier transforms with 
respect to the spatial coordinates. Taking into account that, in this general development, 
the initial conditions are arbitrary functions of the spatial coordinates, the best way to give 
this general solution will be using the convolution product between functions (denoted by 
*) with respect to the spatial coordinates. That is, the properties of the inverse Fourier 
transforms and of the convolution product allow us to write the solution as: 

Hv, = Qi{r], x) * 0i(x) + Q2{v, x) * 0-(x) 

where Qi{r),x) and Q2{r],x) are the Green functions, that is the inverse Fourier transforms 
of the fundamental solutions (f)i{r), s) and 02(^) respectively. Then, we need to calculate a 
few inverse Fourier transforms to obtain the Green functions. 

Firstly, we consider that the fundamental functions (28) and (29) can also be written in 
the following form: 

/ X 3_ / sin kq\ riA3ri — rii) sinkq rif ^ /sm kq\ , 
? / X 9?7i 1 nr^r>\ sin/co 3r7? f sin kq\ nf sinkq 

where g = e (^1 — e = r/rji, k = ■ s and being the operator Dg{f) = f — gdgf. 
This form simplifies the number of inverse Fourier tranforms that we have to obtain. In fact, 
we shall only need a pair of well-known inverse Fourier transforms as we shall see below. 

Let us remind the definition of the inverse Fourier transform of a function with respect 
to the spatial coordinates, that is: 



With this definition, we have that the inverse Fourier transform of the function {sin kg) /k is 
known and has the following general form: 



sin kX 
k 



^ {SD{X-r)H{X) + 6D{X + r)H{-X)} (B3) 



AttX 



-25- 



where 5d{x) represents the Dirac delta distribution and H{X) is the Heaviside unity function. 
On the other hand, we also have the general expression (Guelfand & Chilov 1962): 



Ir/.-A-nl 



2A+nr(^)7r3/2 



(B4) 



where n denotes the dimension of the space where are realized the inverse Fourier transforms 
and r(-) represents the Gamma function. These two expressions will allow us to calculate 
all the inverse Fourier transforms involved in the Green's functions. 

Having a look to the fundamental solutions (Bl) and (B2) we can see that we need the 
inverse Fourier tranforms of functions of the form (sin kg) /k^, which can be obtained from 
(B3) and (B4) using the convolution product in the following way: 



-1 



sin kX 
kP 



sin k\ 
k 



As the convolution product in general is given by the expression: 



f{t)g{x - t)dt 



(B5) 



it results that the two inverse Fourier tranforms needed are expressed in general as: 



1 


sin k\ 






k^ 


47r \ 



H(X -r) + -H(r - X) 
r 



H(-X-r) - -H(X + r) 
r 



H{-X) 



^1 = \r{r^ + 3A^)i/(A - r) + X{}? + Zt')H{t - A)l i/(A)- 

k^ J z47rr L 

- [r(r' + 3A')//(-A - r) - A(A' + 3r')//(A + r)] //(-A)} 

The corresponding operators acting over these expressions will give us that the Green 
functions have the form: 



47re3' 



47re77^ 



9 



vf Soir-g) 



To obtain these expressions we have considered that g is always a positive number (that is, 
H{g) = 1 and H{-g) = 0). 
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Finally, as we have said above, the general solution of the initial value Cauchy problem 
in the real space is given by the convolution product between the Green functions and the 
initial conditions, that is: 

Hv, x) = Qiiv, x) * (j)i{x) + Q2{v, x) * (j)[{x). 

This expression reduces to unidimensional integrals in the case of spherical symmetry. That 
is, if we consider the initial conditions = /i(r) and (/)[{x) = f2{r) as functions depending 

only on the radial coordinate r, then the corresponding convolution product, defined by (B5), 
is written as unidimensional integrals as we show in the following. 

To make clear the expressions we can consider firstly the case when /i(r) ^ and 
f2{r) = 0. In this case the final gravitational potential will be written as: 

'Piv^r) = ^ j q^fi{q)dqH{g-r) + — j qfi{q)[2rq + K{r,q,g)jdqj + 



Vi9i^V Vi) J'^ f,{^r^ + g2_2rgx)dx + ^ J'^ My^r^ + g^ - 2rgx)dx+ 



2 /•! 



4 Ct^B^/'^"^^^^^"^ 



2rf J_i _|_ ^2 _ 2rgx 

where K{r, q, g) — g"^ — r"^ — q^. 

On the other hand, the case in which we have /i(r) = and f2{r) 7^ we will have that 
the final evolution of ^(r^, x) is: 

''^^%{^jy ^'/2(g)i^(r, q, g)dqH{g - r) - ^ qh{q) {K{r, q, gf + ATY)dq+ 

+ - / qK{r,q,g)f2{q){2rq + K{r,q,g))dq\ + ^ f2{Vr' + g' - 2rgx)dx+ 

J\r-g\ ) ^ J-1 

3„2 r rg-r 1 pr+g 

+ ^|2 q^f,(q)dqH{g-r) + - J qf2{q) {2rq + K{r,q, g)) dqj 
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